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Summary 

Algebraic expressions in terms of elementary mathematical 
functions are derived for power absorption and dissipation by 
eddy currents and magnetic hysteresis in ferrite spheres. Skin 
depth is determined by using a variable inner radius in 
descriptive integral equations. 

Numerical results are presented for sphere diameters less 
than one wavelength. A generalized power absorption param- 
eter for both eddy currents and hysteresis is expressed in terms 
of the independent parameters involving wave frequency, 
sphere radius, resistivity, and complex permeability. 

In general, the hysteresis phenomenon has a greater sensi- 
tivity to these independent parameters than do eddy currents 
over the ranges of independent parameters studied herein. 


Introduction 

Damping of electromagnetic waves by iron oxides differs 
distinctly from that by iron (refs. 1 and 2). At low wave 
frequencies large persistent eddy currents can be induced in 
iron because of its high electrical conductivity, a process that 
converts field energy into heat as ohmic losses. At high wave 
frequencies, however, the oscillating currents cannot 
effectively penetrate the surface of iron. 

Ferromagnetic and ferrimagnetic oxides (ref. 3, pp. 24 
and 25), on the other hand, with their low conductivity and 
large skin depth (ref. 4) can be effective in absorbing field 
energy by eddy currents even at high wave frequencies. Their 
low weight densities, approximately 60 percent that of iron 
(ref. 2), make them attractive for flight applications. 

Associated with a harmonic wave motion is a relatively large 
hysteresis effect due to frictional resistance in the crystalline 
structure (Weiss domains) and to radiation from electron spin 
reorientation (ref. 5, p. 16), again resulting in dissipation of 
the applied field energy. These damping mechanisms cause 
a phase shift between the applied magnetic field B A and the 
resulting magnetization M, which enters the analysis by use 
of a complex susceptibility, x — x' — jx" ( re f- 5, p. 100). 
Dissipation of field energy in the study reported herein 
was treated with a complex isotropic permeability 
pt = ^ -jfi" = 1 + x' - jx " in the ferrite spheres. 

Large ratios of wavelength X to sphere radius a are of main 
interest. Scattering is then of the Rayleigh type, which is 
proportional to sphere diameter divided by wavelength to the 


fourth power (ref. 6, p. 457). Therefore diffraction can be 
neglected. 


Analysis 

Power absorbed per sphere due to internal eddy currents 
P ec plus that due to hysteresis P hy is 

+ ri-idV+'-^H-^dV ( 1 ) 

The linear nature of Maxwell’s equations permits P cc and P h> , 
to be treated as separate individual events. Symbols are defined 
in appendix D. 

The objective now is to evaluate these two integrals in terms 
of elementary algebraic functions and the basic parameters of 
the applied wave along with the material properties of the 
ferrites. 

Wave-Particle Interaction 

Consider an applied transverse electromagnetic wave with 
a propagation vector k and mutually orthogonal electric and 
magnetic field vectors and (fig. 1) having a harmonic 
time dependence e jut . The magnetic vector potential relation 
V X A = B used with the applied wave motion of figure 1 
gives 


A A = <p B A r sin 0 (2) 

The large \!a values of interest herein permit neglect of 
Maxwell’s displacement current inside the spheres (ref. 7, 
p. 368). Then with use of an internal magnetic potential A, 
and the laws of Faraday, Ohm, and Ampere, 

V 2 Aj = jo)fir ~ 1 A ; 


With the assumption that A, is independent of azimuthal angle 
<p (as on p. 374 of ref. 7), this equation reduces to 


• - (? S A) 

r 2 dr \ dr J 


(1 -cos 2 #)^ d 2 [(l -cos 2 0)‘^A 


d cos 2 d 


= jo)fiT l A, (3) 


l 


(7) 



Figure L— Transverse wave interacting with a sphere. Wavelength X much 
greater than sphere radius a. 


A 0 = J 1 B a (r + Dr 2 ) sin 0 


Here conductivity is assumed to be negligible, and permeability 
outside the sphere (i a is considered to be real. 

The constants C and D of equations (4) and (7) are next to 
be determined. Solutions must be matched at the surface r — a 
so that (ref. 7, p. 376). 

\ 0 (a)=\ i (a) (8) 

To satisfy conservation of magnetic flux across this boundary; 
B • n must equal zero and thus by Stokes’ theorem 


d(r\ 0 ) 3(rAJ) 

li — = \x 0 — at r = a 

or dr 


(9) 


This equation implies that there is no surface current and that 
the current enters the sphere as a continuous distribution, which 
is the true case for microscopic detail. 

Substituting equations (4) and (7) into both equations (8) 
and (9) gives, in terms of Bessel functions, 


Solving equation (3) by separation of variables gives 


A, = <p ^ B A Cr l/2 / V2 (rV/>) sin 0 


(4) 


where 


a 3 + D = a 3/2 C 


‘ - 1/2 


(a\f]p) 


rjjp 


1/2 


(a'Jjp) 


( 10 ) 


P = 


LOfl 

T 


(5) and 


and 


(2a 3 - D) n = — m ^ 3/2 C 


I- l/2 (a\[jp) 


P " p t ^7p*~ ' ( 6 ) 

By convention, a negative sign is used before j in the definition 
of complex permeability p. The imaginary part \i " must then 
be a positive number for a loss of wave energy. 

The contribution of the applied wave plus the vanishing 
influence of the sphere on the external vector potential with 
radial distance is 


- [WJP + 


ayfjp 


t) 


ha (a'fjp) 

( 11 ) 


Equations (10) and (1 1) can be combined to solve for C and 
D as 


and 


C = 


3na 5/2 \fjp 


{p-Po)asfjp I 1 +jpa 2 ) ~n I U2 {asfjp) 


a'Jjp (2n + p 0 ) I_ l/2 (a\fjp)~ 

[p 0 (l+jpa 2 )+2i. 

*] hn( a 'JJp) 

a'Jjp (p~Po) I-\/2( a 'Jjp) + [poO+jpa 2 )-p] 

h/2 (aJJp) 


( 12 ) 


( 13 ) 
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Thus the equations to determine the potentials A, and \ 0 keep 
the same form as in reference 7, to which they reduce as p" 
goes to zero. 

Eddy Currents 

By using Ohm’s law in Faraday’s equation 

aB 

V X E = rV X i = = -JwV XA| 

dt 

a solution of which, for current i inside a sphere, is 

-M A < 


and 


- W 


(14) 


The power absorbed per unit volume of the sphere by eddy 
currents is 


dP c c 1 . - 

— — = - n* i 
dV 2 


(15) 


c r = cos(2 r dm\fjp) (22) 

in the upper and lower limits of integration, respectively, and 

CC = 9 ir \v 2 1 |p | V2 a b j (ul | p |a 2 ( C a + c a ) 

+ [u 2 g + aVo (a 2 |p | 2 - 2 3m p)] (C a - c a ) - 2 a 
(S a + S a ) (Re“\J^— (S a —S a ) 


u l g + 


2J 


(S a - s a ) (Re [ Up a/| 


+ 2 a p Q 
+ (S a + s a ) 3m l UpyJ r - 




(23) 


where 


dV = r 2 sin 0 dd dr d<p 


(16) 


Substituting equations (4) and (14) into equation (15) and 
integrating the result over 6 from 0 to 7r gives the power 
absorbed per sphere as 


Pec = — BaCC [ hnirtfp) hair^-jP) r dr < 17 ) 
3 T 


if the radial distance to the lower limit of integration d is set 
equal to zero. Varying d allows the amount of power 
absorption versus wave penetration distance into the sphere 
to be determined for a study of “skin depth.” 

The integration is performed in appendix A and the constant 
CC is determined in appendix B. This results in 


= 


u l B z A CC 
3tV|p| (Re p 


(S a - S d ) 3m\^ + (s a - S d ) (R e'JJp 


where U = p - and Uj = pp - 2 p 0 (Re p + p 2 . These 
equations reduce to those of reference 7 as distance d and 
permeability p" go to zero. 

Hysteresis Loss 

Energy of the applied electromagnetic wave is also dissipated 
in the continual attempt of the applied oscillating wave motion 
to alter the Weiss crystalline structure and the electron spin 
orientations inside the ferrite spheres. This is illustrated by 
the lag in B to the changes of H in a hysteresis loop. 

The energy change per unit time P hy per unit volume in a 
hysteresis loop is 

H • — = - jw H • B 

dt 2 

for a harmonic time dependence. For an isotropic permeability 
B = /IH = (p ' + jp") H Thus the wave power loss per unit 
volume due to hysteresis is 


/IP 1 _ 

= - H*H) 


fc a — c a C d — cA <Re p 1 

(18) 

dv L 

\ a d ) \p\ \ 


^ rr tt2 

= - up H 

a and d in 
S r = sinh (2r (Re'JJp ) 

(19) 

2 

-2| „“( |B ' |2+ft|2 ) <24> 

s r = sin (2r 3m V/p) 

(20) 

C r = cosh (2 r (Re'Jjp) 

(21) 

which is consistent with that on page 100 of reference 5. By 
using equation (4) 
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„ A CB, 

B = v x A - = 372 


cos e h a {rsfjp) 

2 sin 6 [ha (rfjp) - hairsfjp) 
0 


Substituting equations (16) and (25) into equation (24) and 
integrating over d<p gives 


P _ pi r r 

r hy ~n~ 


a r T ^o 2 


d J 0 


cos 2 #-f U sin 2 (9 


ha(rjjp) 


X ^3/2 (rsA -“ 7 / 7 ) sin # (id dr 


+ , Ip 1 1 r h/2( rs f]p) ~jp) sin 3 0 dd dr 

4 J o 


-1 Iin(rJjp) h/ 2 (rj-jp) + 

4 Jo L 

x ~jp) hn( r hjp) \ sin 3 0 dr j (26 


Using 


(cos 2 # + 14 sin 2 #) sin 6 dd = 1 


and 


u: 


sin 3 # d6 = 


1 


in equation (26) gives 


^hy “ 


7T60/X 


I/* I 


Bice 


2 ^ 


~^Jjp hn ( r '//p) ha —JP ) '/jp dr 


1 


+ . r ^~]p hai’^/Jp) h/ 2 (r^ —jp) 4jp dr 

i Jd 


(25) 


1 


*/jP A /2 (r^jP) ^3/2 jp ) 


X /i/ 2 (rV-y» / 3 / 2 (rv(^)J dr 


(27) 


The integrations over rare shown in appendix C. The resulting 
expression (Cl 8 ) can be written as 


B hy = 


0)JJL 

H 1 

1 

+ - 
3 




« 2 + ff 2 t \ C r 4- c r 


S\ 7V 


aS r + I3s r C r — c r S, _ 
(Vjpjr ) 3 a 0 


w 2 


( 4 0 2 a 2 \ 

+ (- + — - y ) a Shi(2r <Res//p) 


/ 4 a 2 0 2 \ _ 

“ (3 + J ~ j) ^ Si(2r 3 m J]p) 


(28) 


/-=</ 


where use was made of definitions (19) to ( 22 ) and 

2 ( ReV/p 

a = 

\P\ 


0 = 


2 tj m'jjp 

Ip I 


(29) 


(30) 


The hyperbolic sine and sine integrals can be expressed as 


Shi(X) = X + — — + 


+ 


and 


Si(A) = X - 


3x3! 5x5! 7x7! 


+ 


(31) 


3x3/ 5x5/ 7 x 7/ 


(32) 


(ref. 8 , eqs. 672.11 and 431.11). 

Generalized Parameters 

The dimensionless power parameter tpI (w Vfii) for both 
eddy currents and hysteresis can be written in terms of two 
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independent dimensionless parameters: a normalized 

permeability 


P -w 
= 

Mo 


and a combination of the remaining physical variables 


Results and Discussion 

For physical clarity it is advantageous to first relate results 
directly to the basic independent variables r, u, m> an ^ r as 
well as to a in figures 2 to 4. The ranges of generalized 
parameters and a/\ are indicated in the legends of these 
figures. The plotting variables in figures 5 and 6 are 
dimensionless power absorption parameters that are functions 



5 






Figure 3.— Effect of varying wave frequency / with radius of sphere 
Q “ 10 m; resistivity 7 = 10 STm; and normalized permeability 
= 100 —j 10. The curves reach asymptotes at / < 10 9 for hysteresis 
and/ < 3x 10 “ for eddy currents. Parameters a/\ — 3.3 X 10~ 13 /and 
t = 8t 2 X 10“ 16 /. 



Figure 4.— Effect of varying radius of sphere a from 2x10 4 to 
256 x 10 m with normalizing radius a a = 10 4 m; wave frequency 
/= 10 Hz; resistivity 7 = 10 ITm; and normalized permeability u/(x () 
= l00-jt0. Para meter a/\ ranges from 6.7X10 -4 to 8.6xl0~ 2 and 
parameter 0t u from 3.2 xlO -5 to 0.5. 


of real and imaginary parts of permeability for fixed values 
of a 0 . 

Skin Depth 

For good absorption of wave energy the electromagnetic 
waves must penetrate the spheres to a sizable depth. The 
fraction of total power absorbed by a sphere as a function of 
the lower radial integration limit r = d in equations (18) and 
(28) is used as an indication of skin depth in figures 2 to 4. 

Figure 2 shows the effects of resistivity r ranging from a 
typical value for iron of 10 ~ 7 Q-m to values that lie on the 
asymptotic curve. This curve was reached as r increased 
beyond 10 _4 il-m for eddy currents and beyond 1 for 
hysteresis. Held constant were frequency at 10 7 Hz, 
permeability ratio at lOO-ylO, and radius at 10 “ 4 m. In 
general the wave penetration depth was considerably greater 
for hysteresis than for eddy currents. 

Decreasing frequency /= u/2it gave the same trends in 
figure 3 as increasing t, as expected from the basic parameter 
a 2 0)(i o /T. Here r held constant at TO Q-m. The ratio a/\ 
varied from approximately 1 to less than 10 “ 3 in the curves 
shown for hysteresis. For eddy currents, however, a/X< I 
yielded curves that lie on their asymptote. Curves for a > \ 
looked reasonable but violated the neglect-of-displacement- 
current assumption and therefore are not shown. 

The curves of various size spheres in figure 4 are for radii 
less than 0. 1 wavelength and therefore are well within the 
validity of the neglect-of-displacement-current assumption. 
Radius was varied by multiples of 2 from 10 ~ 4 to 
256 X 10 “ 4 m. Again t and w were held constant at 10 and 
2ir x 10 9 , respectively. The eddy current curves, shown 
dashed, are asymptotic at a<3.2xl0~ 3 ; the hysteresis 
curves are asymptotic at a <2 X 10 _4 . 

Power Absorptions and *Vy 

The same parameter, r/a 5 u 2 B % factors out of both the 
hysteresis equation (28) and the eddy current equation (18). 
Typical results, shown in figures 5 and 6 for these cases, are, 
however, quite different. 

The eddy current power parameter was insensitive to d 0 
over the wide range of 10 ~ 7 < d 0 < 10 ~ 2 as shown in 
figures 5(a) and (b), but gradually decreased with increasing 
d a beyond 10 ” 2 , as shown in figures 5(c) and (d). These 
three-dimensional plots are symmetric about a diagonal curve 
(shown dashed) where fi' — . This curve also represents 

a local minimum of the eddy current power parameter at low 
f. i' and \i n . 

The eddy current power parameter varied over less than 3 
decades, but the hysteresis power parameter varied from 12 
to 7 decades, depending on (t^; as /z' and \x ff were varied 
over 3 decades in figures 5 and 6. 
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Figure 6. — Concluded. 
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Real and imaginary permeability played similar roles in the 
eddy current formulation, since they entered P cc only through 
p. Starting from the dashed diagonal line of figure 5 the 
parameter TP eQ /a 5 o) 2 B A increased to the same peak values 
with increasing p' or p" . 

Power absorption by hysteresis, however, was directly 
proportional to a multiplier p" in addition to its dependence 
on p. This offset the symmetry of dependence on p f and p" 
present through the parameter p in equations (18) and (28) 
and made absorption much more sensitive to hysteresis than 
to eddy currents. As a result an upper elevation (shown dashed) 
of the three-dimensional plots of TP hy ia 5 u 2 B A formed 
close to the diagonal p" /p 0 = p'lp 0 in the range 

io _7 ^a 0 <io _I . 

As a practical example consider a resistance r of 
10 4 12-m, a B a of 1 gauss (10 “ 4 tesla), and a frequency a? of 
2tt X 10 10 rad/sec. Then the peak power per sphere absorbed 
by hysteresis varies from about 0.2 to 14 W as the sphere 
radius a varies from 0.1 mm to 1 cm in figures 6(a) to (e). 
On the other hand the maximum power per sphere absorbed 
by eddy currents varies from about 1.4 X 10~ 10 to 1.4 W for 
the same r, B A , and range of a in figures 5(a) and (b). 

For this same example let p ” approach zero in figure 5 while 
holding p ' equal to p a . The power per unit volume by eddy 
currents then reduces to 0.2 W/cm 3 for the same r, B A , and 
o) values. This power absorption is in close agreement with 
results obtained from the sphere equation on page 135 of 
reference 5. 

Concluding Remarks 

Expressions were obtained, in terms of simple algebraic 
functions, for power absorption of electromagnetic waves in 
ferrite spheres by eddy currents and hysteresis. Skin depth was 
especially large for hysteresis and was significant for eddy 
currents, indicating good absorption for the small spheres of 
special interest herein. 


Generalized parameters were found that reduced the five 
independent variables wave frequency o?/27r, sphere radius a, 
resistivity r, and real p f and imaginary p" permeability 
divided by free-space permeability p 0 to three independent 
generalized parameters d 0 = a 2 up 0 /T, p'/p Q , and p” !p a 
and two power loss parameters TP ec ia 5 io 2 B A and rP hy fa 5 o: 2 B A 
for small spheres. Working curves were presented for 
obtaining power losses from input to the independent 
parameters. 

The eddy current power parameter was insensitive to Q 0 
over wide ranges, 10 ~ 7 < Q 0 < 10 ~ 2 , and gradually 
decreased with increasing <2 0 beyond 10 “ 2 . Real and 
imaginary permeability played similar roles in eddy current 
formulation, since they entered P ec only through complex 
p~p' — jp" - Power absorption by hysteresis, however, was 
directly proportional to p" in addition to its dependence on 
complex p. This offset the symmetry of dependence on p' 
and p " and made absorption much more sensitive to hysteresis 
than to eddy currents. 

From results obtained by using the equations derived, it 
appears that ferrite spheres provide a good means of damping 
extraneous electromagnetic waves and isolating electronic 
circuitry, especially in flight applications where weight is a 
consideration. 

The components of the induced magnetic field, B r and B e , 
diminished rapidly with radial distance as r -3 and r~ 2 , 
respectively. Except for closely packed spheres, there should 
then be little interaction between small spheres suspended in 
an insulator. The intentional degradation of an electromagnetic 
wave by its penetration into such a medium should then be 
closely proportional to the distance of wave travel and the 
number density of spheres. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, June 13, 1989 
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Appendix A 

Evaluation of Integral in Equation (17) for 
Power Loss by Eddy Currents 


The equation of the modified Bessel function l m (kr) can be 
written as 


1 d 
~rJr 


dl m (kr) 

dr 


= U 2 + — ) Ukr) (AD 

r 


and its conjugate as 

df m (kr) 


1 ^ 

r dr 


dr 


k + ~2 I l m(kr) 
r 


where 


k~'Jjp={ \ 
k = 4^jp = (i 


(A2) 


(A3) 


Multiplying equation (Al) by I m (kr) and equation (A2) by 
I m (kr ) and forming the integral of the second of these 
equations subtracted from the first gives 


- d 
L(kr)- 

dr 


dl„,(kr) 

dr 


-Im(kr) - 

dr 


dl m (kr) 


dr 


dr 


= j(p+P ) l I m (kr) I m {kr)r dr 

' d 


Using integration by parts on the left side yields 


, , T , dl m (kr) dl m (kr) 

I m (kr)r : l m (kr)r 


dr 


dr 


j(p+p) \ I,„(kr)I m (kr)r dr 
*d 


Finally 

jip+P ) j I,„(kr)I m (kr)r dr 


r I m (kr)kl^(kr) - I m {kr)kl,’ n (kr) 


(A4) 


where the prime marks denote differentiation with respect to 
the arguments of the Bessel functions. 

Substituting equation (A4) into equation (17) with m - 3/2 
gives 


p = 

1 ec 


nrBjCC 
3 ]t(p+p) 


kl y2 (kr)li /2 (kr) 

- kly 2 {kr)ly 2 (kr) 


(A5) 


' a 

r dl m (kr) dl m (kr) 

dl m (kr) dl m (kr) 

J d 

dr dr 

r 

dr dr 


r dr = 0 


Using m — 3/2 in 


m 


KM) =fm- lU) --ImU) 

Z 


(A6) 


of reference 9 (p. 79, eq. (3)) to eliminate the derivatives in 
equation (A5) yields 


p = 

1 ec 


3 jr(p+p) 


r \ kI i/2 (kr)I V2 (kr) 

— kK/ 2 (kr)I \/ 2 (kr) 


Then using m = 1/2 in equation (A6) and reference 9 (p. 79, 
eq. (1» gives 


r—d 
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P' c = 


7tu rB 2 A CC 


3 hip+p) 

\kl_i /2 (kr)l m (kr) - kl. m (kr)l u2 (kr) 
j(i P+P) 


Ip I 


hn{kr)l xn (kr) 


which can be written as 


p = 

1 ec 


2tt o) 2 g^CC 
3r(/>+p) 


3m j^rW_ j /2 (^r)/[/2(fcr) 

(Re p 


Ip 


hn&rWmite) 


(A7) 


Utilizing the relations 


r x 2 cosh z sinh z 

I-\a{z) I\ /2 (z) = — 

7T k ! 


^2(ReZ _ ^ __ ^ — 2/£W 


2t|z | 27T fz I 

= —j^-j ^sinh(2 (Re z) + j sin(2 3*n z)J 

(A8) 


and 


, . . . 2 sinh z sinh ; 

hnUVuiiz) = — 

7T |z| 


„2(R tf + £ ,-29i t ; g 2j3~z , - 2/3, 


2x|z 


27t|z ! 


= — — |cosh(2 (Re z) — cos(2 3m z)J 


(A9) 


in equation (A7) with z = kr = \fjp r and z — kr gives 

2u 2 B 2 a CC 


P = 

1 ec 


3 r(p+p) 



JP 


\P\ 


sinh(2r (ReV//? ) 


1 (Rep 


1 3/2 


+ y sin(2r 3mV^)Jj ^ ^ 

X j^cosh(2r (Re'Ijp) - cos(2 r £JmVj^)J 
co 2 ^CC 


•%/j — j y 


3 T (Re p 

X ^sinh(2r (ReV/p ) 4- j sin(2r Sm'Jjp)^ 


(Re p r n- ■ 

— p /2 [cosh(2r (Re'Jjpy cos(2 r £JmV//?) 


Finally 


/> = 


a) 2 fljcc 


(Sm'Jjp) sinh(2r (Re'Jjp ) 


3tV[p| (Re p 
+ ((ReV^) sin(2r 3mVjp) - ^-E. 


»• IP I 


X j^cosh(2r (Re'JJp ) — cos(2r 3mV^)j 


(A 10) 


For d — 0 and fx " — 0 (therefore p = (Re/?) , equation (A 10) 
reduces to that of reference 7 (p. 378): 


p = 

* ec 


u^CC 

~3rp 3/2 a 


- (5 + J) - (C-c) 


where w = aV2p, S = sinh p, s = sin fx, C ~ cosh 
c = cos and, by using equation (A3), EfmV/p = (ReV^ 

= >fp/2. 
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Appendix B 

Determination of Constant CC Used in Equations (18) and (28) 


The scalar product of constant C in equation (12) times its 
complex conjugate can be expressed as 


CC 


9 fifi 4pp a 5 
AA + BB T ab + AB 


(Bl) 


where the constants in the denominator are defined by 
AA = jV/I - 2fi 0 Re \i + 

X \p\a 2 l- xn {'ljp a)I_ l/2 (/~jp a) (B2) 


BB =[p- Ho ( 1 + jpa 2 )\ [p - p 0 ( 1 - jpa 2 ) 

X /i /2 (> Jjp a)I u2 ('t-JP a) (B3) 

AB = (n - p 0 )a 'JJp [/* 0 (1 -jpa 2 ) ~ p] 

X I-vii'fjP a) f m (^]p a) 


and 

AB = (£ - p„ ) a V rz # [/i„(l +jpa 2 ) - m] 

x /_ | /2 ( S-Jpa ) / 1/2 (\f]p a) (B5) 


I-vii'JjP «) ^ — 1/2 «) 


1 

ir/~\p\ a 

X j^cosh(2 Re/jp a) + cos(2 3m-J]p a) 

(B6) 


Since 

a ) A/2 (nA -jp a) 


= 1-inS'T-jP a ) T \ M'fip a ) 

it follows from equation (A8) that 


I-irti'JJp a) h/jt'f-Jp a ) 

1 

ira/lp | 

X [sinh(2tf Re '/jp) -j sin(2 a 5 m -/jp) 

(B7) 


The product of Bessel functions in equation (B2) can be written 
by using reference 9 as 


I-\nS'ljpa) I-\n('f-jP a ) =' 


-k VH 


cos 


ih(V/p a) 


X cosh {'J—jp a) 


Using equations (A9) and (B7) in equations (B3) and (B4), 
respectively, along with equations (B6) and (A8) in equations 
(B2) and (B5) reduces the Bessel function expressions in 
equations (B2) to (B5) to simple hyperbolic and trigonometric 
functions defined in equations (19) to (23). This results in the 
expressions 


1 


(e^ pa + e~' Tpa ) 


AA = a-Jty j U: 


C —c 

i '-'a 


(B8) 


2ir\f\p-\a 

X ( e ' p ^ a + e ~'^jp a ) 


BB = [ U 2 g + fi 2 ( a 4 \p | 2 - 2a 2 3m />)] ( B9 > 


and therefore 


and 
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N 


■Uj [(5 a + s a ) (Re v|- (S a - sj 3m V|] 

+ [ (S a -s a ) (Re (Vpyfa + ($, + s a ) 3m (t/pV^)]) (BIO) 


Substituting equations (B8) to (BIO) into equation (Bl) gives 
equation (23). 
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Appendix C 

Evaluation of Integrals in Equation (27) 
for Power Loss by Hysteresis 


Equation (27) can be written as 


1 T - - z 1 , - -z, 

3, = — e- + e z (e- — e *) 

2r L 2 


ImI 2 ^i 


noifi " «;,CC 


%—-6 l + - (3 2 - 3 3 - 3 4 ) 


(Cl) 


x[ 


e z + e z — - (e z — e ") 

z 


dz 

zk\ 


Here 


= i“ ^3/2(2) ^3/2(2) dz 


and by using reference 10 (p. 443) 

d,=- I — ( cosh z-- sinh ) fcosh z~z. sinh z) dz 

xJz|z| V z M * ' 


In like manner 


2 T 2 z 

S 2 = - zI m (z)r W 2 (z) dz = - — sinh z sinh z dz 

7 T ft . K 


3 3 = U/ 2 (z) ^3/2 (?) dz 


= - \ — sinh z fcosh z - — sinh z) dz 
Trjlzl . V Z > 


and 


2 1 _ / 1 

— — sinh z , _ 

x |z| v z 


<J 4 — £J 3 = - 1 — sinh z (cosh z sinh z) dz 


where 


1 

27T 


'-ir + Hr 


,-iy + ( 1+ iu- 


dz 

zlzl 


2 ir J L\ zz z z) 


,2<R«z 


+ (1 + — + - + 4r) e~ 2<RcZ 
zz z z 

+ , 1 - — - 1 + — ) e 2j3mZ 
ZZ z z, 


1 1 1 

+ [ 1 - — + — -) e 

zz z z 


-2 j9" 


dz_ 

z|z ! 


■i! 


ZZ + l I e 7(R<z + 

I 

ZZ 


+ 


zz 


— 


( 1 ^ + r ^|.£±i(e^-r 2 * i! 


zz 


+ i ( e 2 -' 3 "* - e _2 ^ w 
zz 


z|z| 




—i cosh(2 (Re z) + cos(2 3m z) 

ZZ zz 


z = r 'JJp and z = rV —jp 


(C2) 


Next, replace the hyperbolic functions by their exponential 
forms, rearrange, and reduce as in the following example for 
3,. 


s inh(2 (Re z) 


ZZ 


2 9m Z . , 

— - — sm(2 9m z) 


zz 


dz 

z|z| 


(C3) 
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\ 


In like manner 
1 


n = - I cosh(2 (Re z) - cos(2 <Jm z) 
7 r 




(C4) 


and 




sinh(2 (Re z) —j sin(2 3m z) 


cosh(2 (Re z) - cos(2 3m z) 


zz 


dz (C5) 


Since 3 4 is the complex conjugate of 3 3 , it follows that 

£J 3 + 3 4 = 2 (Re 3 3 (C6) 

Using equations (C2) and (C5) in equation (C6) gives 

3 3 + 3 4 = j sinh(2 (Re z) 

-3msjjp S in(2 3 m z)-^ii# 

kl 

x ^ cosh (2 (Re Z) - cos(2 3m Z)}} — (C7) 


Substituting equations (C3) to (C7) into equation (Cl) and 
collecting terms of like integrals gives 


\dl 2 P hy ( a 3 , (3 2 a 

= ( a H 


CC 


6 6 
cosh x 


, I n P 3 

dx+ [ j8 — — H - 

6 6 


cosh jc 


■I 

-It 1 **-'! 

-'iT-i! 


I , ■, 1 sinh x . 

X l cosy dy - or' \ — — dx + (3 3 


dx 


cosh x dx — ■ 


(C9) 


3/3 


17 ^ 1 - 


JC 

sinh x 


dx 


(C8) 


where 

x = 2 (Re z 
y = 2 z 

2 6ie\fjp 
a = 7 = — 

_ 2 3mV/> 

VR 

The following expressions enable the integration of equation 
(C8). They were obtained by integrating by parts along with 
using expressions from pages 89, 96, 148, and 149 of 
reference 8. 


I 


cosh x cosh x 

— dx = h Shi x 


(CIO) 


1 


cosh jc 


dx = — 


cosh jc 

sinh jc 

cosh x 

3x 3 

6x 2 

6x 

-dy= - 

cos y 

f sin y 

y 

J y 


4- - Shi x (Cl 1) 


(Cl 2) 


cos y cos y sin y cos y 1 

„4-^=--rx + — f + — ^ + -siy 


f cosy 

JT" 

1 


3y 6 y L 6 y 6 


(Cl 3) 


sinh x sinh x cosh x 1 

— -dx=~~— 1 — + - Shi x (Cl 4) 

x i 2x z 2x 2 


It’— 


sin y cos y 1 

— —z I — Si v 

2y 2 2 y 2 


(Cl 5) 


where 


Shi jc - 


and 




sinhx jc 3 x 5 x 1 

dx = x + -4 1- + 

x 3x3! 5x5! 7x7! 


• (Cl 6) 


-if 


c- . Sl °y , J 3 >’ 5 y 7 

— + + (Cl 7) 
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These expressions then enable the integration of equation (C8) 
to the following form: 


iM 2 ! Phy _ 


a 3 a/3' \ cosh x 

— H a r 

6 6 ) x 



cos y 

y 


a 3 sinh x /3 3 sin j 

+ T jc 2 + 3 >’ 2 


a 3 cosh A' j 0 3 cos y 

I x 3 + T 


1 . 1 . (4 B 2 <x 

3a 3/3 V 3 6 





.It >3’2 


i 


(Cl 8) 


where a 2 and y 2 are the real and imaginary parts of z (see eq. 
(C9)) at the outer radius a and A| and y { are at an arbitrary 
inner radius of interest. 
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Appendix D 
Symbols* 


a a 0 n„ 

Q a a 2 p u 

A magnetic vector potential 

a radius of sphere 

B magnetic flux density 

B magnetic field vector 

C constant defined by eq. (12) 

C r constant defined by eq. (21) 

c r constant defined by eq. (22) 

D constant defined by eq. (13) 

3D constant defined by eq. (34) 

d radial distance to lower limit of intergation inside 
sphere 

E electric field vector 

f wave frequency, co/27r 

H magnetic inductance 

3 integral in appendix C 

£fm imaginary part of a complex quantity 

l m modified spherical Bessel function where m is an 
odd integer plus 1/2 

i current 

j 

k yfjp 

M magnetization 

n normal to boundary 

P power 

p parameter c op It 

p D parameter top 0 lr 

(Re real part of a complex quantity 

r radial distance from center of sphere 

S r defined by eq. (19) 

s r defined by eq. (20) 

t time 

U p-p o 

u, I 1 


* All terms are in international units unless stated otherwise. 


Ug M~ — 2 Po (Re P + pi 
u]„ \p„\ 2 + (1-2 (Re p„) 

V volume 

X dummy variable in eqs. (31) and (32) 

x 2 (Re z 

y 2 z 

Z rfjp 

a defined by eq. (29) 

defined by eq. (30) 

9 polar angle 

k wave propagation vector 

X wavelength 

p permeability inside sphere, — jjx" 

p n normalized permeability, p/\i a 

Po permeability outside sphere (vacuum), 47 tX 10 ~ 7 
henry/meter 

t resistivity 

v? azimuthal angle 

X susceptibility 

co angular velocity 

Subscripts: 

A applied wave 

a at a radius equal to a 

d at a radius equal to d 

ec eddy current 

hy hysteresis 

i inside sphere; internal 

o outside sphere 

r radial direction 

6 direction of polar angle 

Superscripts: 

real part of complex quantity 
imaginary part of complex quantity 
complex conjugate 
unit vector 
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